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Abstract 

A statistical meclianical approach to granular material is proposed. Using 
lattice models from standard Statistical Mechanics and results from a mean 
field replica approach we find a jamming transition in granular media closely 
related to the glass transition in super-cooled liquids. These models reproduce 
the logarithmic relaxation in granular compaction and reversible-irreversible 
lines, in agreement with experimental data. The models also exhibit aging 
effects and breakdown of the usual fluctuation dissipation relation. It is shown 
that the glass transition may be responsible for the logarithmic relaxation 
and may be related to the cooperative effects underlying many phenomena of 
granular materials such as the Reynolds transition. 



I. INTRODUCTION 

Despite their importance for industrial applications, non-thermal disordered systems 
as granular media have just recently begun to be systematically studied in the physics 
community In particular, concepts from Statistical Mechanics seems to be successful to 
describe these systems as suggested in his pioneering work by Edwards [0. Actually, granular 
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media are composed by a large number of single grains and, as much as systems of standard 
Statistical Mechanics, each of their macroscopic states corresponds to a huge number of 
microstates. Furthermore, they show very general reproducible macroscopic behaviours 
whose general properties are not material specific and which are statistically characterised 
by very few control parameters, such as their density, load or amplitude of external drive [|I| . 
Granular media are "non-thermal" systems since thermal energy plays no role with respect, 
for instance, to gravitational energies involved in grain displacements, however, thermal 
motion can be replaced by agitation induced by shaking or other driving and this allows the 
system to explore its space of configurations. 

A "tapping" dynamics has several control parameters. A very important one, related to 
the amplitude of vibrations of the external driving, has appeared to be the adimensional 
ratio Tex = a/g ^j^, where a is the shake peak acceleration and g the gravity acceleration 
constant. It has been suggested that Tex plays a role very similar to "temperature" in thermal 
systems. For simplicity, below we suppose the characteristic frequency, u, of external driving 
is fixed, and actually in typical experiments as those in Ref.s we are interested in, it 
plays a minor role with respect to T^x (see [|I|). Material parameters such as restitution, 
friction coefficients, presence of moisture, etc., are also to be usually considered, but they 
generally do not affect the overall behaviours. 

Experimentally it is found that at low T^x relaxation processes as density compaction 
under tapping are logarithmically slow. An essential ingredient in the dynamics of dense 
granular media is, instead, the presence of "frustrating" mechanisms on the motion of grains 
due to steric hindrance and hard core repulsion. We have introduced Frustrated Lattice 
Gases to describe the dynamics of granular materials in the regime of high densities or not 
too strong shaking . They interestingly allow to interpret, in a single framework, several 
different properties of granular media ranging from logarithmic compaction, density fluctu- 
ations, "irreversible-reversible" cycles, aging, breakdown of fluctuation dissipation relation, 
segregation, to avalanches effects, the Reynolds transition and several others P,|^,|TB|,IT3,P5 



In particular they predict the existence of a definite "jamming" transition in granular media. 



in strict correspondence with the glass transition of glass formers and spin glasses |ll|,|8|-[l0| . 
Furthermore, they are new important applications of the statistical mechanics to powders 
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IL FRUSTRATED LATTICE GAS MODELS FOR GRANULAR MEDIA 

The very simple schematic models for gently shaken granular media we introduced are 
based on a drastic reduction of the degrees of freedom of the systems to those we suppose 
are essential. The models consist of a system of elongated grains which move on a lattice. 
Grains are subject to gravity and, eventually, to shaking, which is simulated with a driven 
diffusion like Monte Carlo dynamics. The crucial ingredient of the models is the presence 
of geometric constraints on the motion of the grains, thus, they have been called Frustrated 
Lattice Gases Particles in our models occupy the sites of, say, a square (or cubic) 

lattice (see Fig.|^). They have also an internal degree of freedom Si = ±1 corresponding to 
their possible "orientations" along the two lattice axes. Two nearest neighbour sites can be 
both occupied only if the particles do not overlap. 

Each "experimental" tap can be divided in two parts: one where the grains average 
kinetic energy is finite and a second where it goes to zero. Thus, in our models grains 
undergo a schematic driven diffusive Monte Carlo (MC) dynamics: in absence of vibrations 
they are subject only to gravity and they can move just downwards always fulfilling the 
non overlap condition; the presence of vibration is introduced by allowing the particles to 
diffuse with a probability pup to move upwards and a probability pdown = 1 — Pup to move 
downwards. The quantity — ln(xo)/2, with xq = Pup/Pdown, as we will see, plays the role 
of an effective temperature and can be related to the experimental tap vibrations intensity, 
r 

ex- 
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A. Hamiltonian description 



The general model introduced above can be described in terms of a standard lattice gas 
of Statistical Mechanics [^0. The system Hamiltonian must have an hard core repulsion 
term (J — >■ oo): 

Hhc = J^fij{SuSj)ninj ■ (1) 

{ii) 

Here rij = 0, 1 are occupancy variables describing the positions of grains, Si = ±1 are "spin" 
variables associated to the orientations of the particles, J represents the infinite repulsion 
felt by the particles when they have the wrong orientations. The hard core repulsion function 
fij{Si,Sj) is or 1 depending whether the configuration Si,Sj is right (allowed) or wrong 
(not allowed), see Fig.|l[ 

The choice of fij{Si, Sj) depends on the particular model. In particular, here we consider 
two models: the Tetris and the Ising Frustrated Lattice Gas (IFLG). The Tetris model is 
made of elongated particles (see Fig. |l|), which may point in two (orthogonal) directions 
coinciding with the two lattice bond orientations. In this case fij{Si,Sj) is given by 
fij'^^^'^i.Si, Sj) = l/2{SiSj — eij{Si + Sj) + 1) here eij = +1 for bonds along one direction of 
the lattice and e^j = — 1 for bonds on the other. In order to have a non trivial behaviour the 
dynamics of the Tetris has a crucial purely kinetic constraint particles can flip their "spin" 
only if three of their own neighbours are empty. 

A real granular system may contain more disorder due to the presence of a wider grain 
shape distribution or to the absence of a regular underlying lattice. Typically, each grain 
moves in the disordered environment generated by its neighbours. In order to schematically 
consider these effects within the above context, an other kind of "frustrated lattice gas" was 
introduced, made of grains moving in a lattice with quenched geometric disorder. Such a 
model, the Ising Frustrated Lattice Gas (IFLG), has the following hard core repulsion func- 
tion, fij{Si,Sj) 0: fJj^^'^{Si,Sj) = l/2{eijSiSj — 1) where eij = ±1 are quenched random 
interactions associated to the edges of the lattice, describing the fact that particles must sat- 
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isfy the geometric constraint of the environment considered as "practically" quenched (see 
Fig.|l]). The IFLG shows a non trivial dynamics without the necessity to introduce kinetic 
constraints. 

The phase diagram of the Tetris Hamiltonian corresponds to the usual antiferromagnetic 
Ising model with dilution. The Hamiltonian of the IFLG exhibits richer behaviours. In 
the limit where all sites are occupied (n^ = 1 Vi), it becomes equal to the usual ±J Ising 
Edwards- Anderson Spin Glass [|10|- In the limit J oo (which we consider below), a version 
of Site Frustrated Percolation is recovered fT^jl^ . 



The other important contribution to the full Hamiltonian of a granular pack we consider 
below must be gravitational energy: Ti. = Hue + Hq, where Hq = gY^ifT-iUii and g is the 
gravity constant and Ui is the height of particle i (grains mass and lattice spacing are set to 
unity). The temperature, T, of the present Hamiltonian system (with J = oo) is related to 
the ratio xq = Pup/Pdown via the following relation: e"^^/"^ = xq. It is useful to define the 

— 1/2 

adimensional quantity F = 1/ ln(a;o ) = T/g, which we assume plays the same role as the 
amplitude of the vibrations in real granular matter, that is F is a smooth function of Fe^. 



It is easy to show that if the system reaches equilibrium, the temperature T is: 

where Q is the number of configurations corresponding to the gravitational energy E. Note 
that for samples with constant particle density, E is proportional to the volume and thus 
1/T coincides with Edwards' compactivity [0 apart from a proportionality constant. 

III. THE DYNAMICS OF COMPACTION 

To describe experimental observations about grain density relaxation under a sequence 
of taps a logarithmic law was proposed in Ref. 0: 

p(t„) = Poo - Apoo/[l + B \n{tjn + 1)] (3) 



5 



This law has proved to be satisfied very well by relaxation data in the IFLG model 
which can be excellently rescaled with experimental data. In Fig.^ MC compaction data for 
four different amplitudes, F, as well as the experimental data for three different amplitudes, 
Fex, are collapsed on a single curve using eq.(^. The agreement is very satisfactory (details 
about these data are in Ref. ). 

Interestingly, the results from Tetris are similar [0, but the asymptotic density poo is 
numerically indistinguishable from 1, thus almost independent on xq, a fact in contrast with 
both IFLG and experimental results 0. 



IV. GLASSY BEHAVIOURS OF GRANULAR MEDIA 

We have seen that compaction shows extremely long relaxation times. Actually we show 
that granular systems are typically in off equilibrium configurations and that they may un- 
dergo a "jamming" transition. In particular, in the same framework of the above models, we 
show how experimental results about "memory" effects (the so called "irreversible-reversible" 
cycles) outline the existence of a dynamical glassy transition, which can be defined in a sim- 
ilar way to the glassy transition in real glass formers. 



It has been shown in mean field approximation and numerically for finite dimensional 
systems [jl6|, that the IFLG, in absence of gravity, exhibits a spin glass (SG) transition at 
high density (or low temperature) similar to the one found in the p-spin model. J.J. Arenzon 



has extended the mean field solution of the IFLG in presence of gravity showing that 
at low F the system is frozen in a SG-like phase, but at higher F it separates in a frozen SG 
phase at the bottom and a fiuid phase on top. Above a critical F the system is entirely fiuid 



A. "Irreversible-reversible" cycles and the dynamical glass transition 



Tapping experiments typically show "irreversible-reversible" cycles (seeFig.^: during 
a sequence of taps, if the system is successively shaken at increasing vibration amplitudes, 
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its bulk density typically grows, as in compaction, and then, after a characteristic Tex starts 
decreasing. However, if the amplitude of shaking is decreased back to zero, the density 
generally doesn't follow the same path since it keeps growing. In this sense these observations 
outlines the existence of "memory effects" . 

In analogy to real experiments, cycles of taps were performed in the IFLG in which the 
vibration amplitude F was varied in sequences of increments 7 = AF/r, with constant taps 
duration r. After each tap the static bulk density of the system p(F„) {n is the n-th tap 
number) was then measured. The data are qualitatively very similar to those reported in 
real experiments on dry granular packs H. Furthermore, they allow to define a "jamming" 



transition point Fg(7) in analogy to the glass transition in glass formers [T^, as explained 
in Fig.|[ 



B. Two-times correlation functions 



Granular media in the above dynamical situations show "aging" too |]T3[. To see this, as 
much as in glassy systems, it is useful to consider two time correlation functions as {t > t') : 
C{t,t') = [(p(t)p(t')) - (pit)) {pit'))]/ [{pit' f) - {pit')f] where p(t) is the bulk density of the 
system at time t. 

The fit, at low F, of the two time correlation function, C(t,t'), (on five decades in MC 
time), reveals the following interesting approximate scaling form: 

c(tt')-(i-c Mt' + ts)M ^^ .4. 

where r (which now is not the "tap duration" as before), tg and Coo are fit parameters. Very 
interesting is the fact that the above behaviour is found in both the discussed models (Tetris 
and IFLG) [|1^. The data for the two models, for several values of F, rescaled on a single 
universal master function, are plotted in Fig. ^ It is interesting the existence of such a 
scaling behaviours in off equilibrium dynamics of apparently different systems pi|. 
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C. Response functions and the "fluctuation-dissipation" relation 



In Ref. |n] it was argued that in granular media it is possible to formulate a link between 
response and fluctuations, i.e., the analogous of a "fluctuation dissipation theorem" (FDT) 
(see [0), where the amplitude of external vibrations plays the role of the usual "tempera- 
ture" . In typical situations the FDT coincides neither with its usual version at equilibrium 
nor with the extensions valid in off equilibrium thermal systems in the so called "small en- 
tropy production" limit |2T|. The origin of universalities in off equilibrium dynamics in the 
above heterogeneous classes of materials is still unanswered. 

Our systems are "shaken" at a given amplitude xq and the average grains height, /io(^) = 
{H{t)) (with H(t) = Y^iViit)/^ )) recorded as long as the "mean square displacement" 
B(t,t') = {[H(t) —Hit'yf'). To measure the response function, the average height, hi(t,t^), 
was also recorded in an identical copy of the system (a "replica") perturbed by a small 
increase of the shaking amplitude, after a time t^. 

The difference in heights between the perturbed and unperturbed systems Ah(t,tiu) = 
hi{t,tw) — hQ{t) is by definition the integrated response. FDT or its generalisations concern 
the relation between Ah{t, t^) and the displacement, B{t, t^). In analogy to thermal systems 



0, in Ref. was proposed that: 

A/i(t,t^)~yA(r-i)B(t,t^) (5) 

Eq.(^ states that, after transients, in a granular system the measure of the height variation 
after a change in shaking amplitude, should be proportional to the "displacement" recorded 
during an unperturbed run. In the simplest situations, as equilibrium thermal systems, 
the proportionality factor, X, is a constant, but more generally, X is function of tw and 
t themselves. Interestingly, in glassy systems the quantity gT/X has the meaning of an 
"effective temperature" p2| , p3| . 

In the present model, eq. (^ seems to be approximately valid, also if in typical off- 
equilibrium situations X slowly depends on F and tu,. This is shown, for high xq, in the 
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top panels of Fig.^ in the long time regime X is equal to 1, showing that, in the high 
"temperature" and low density region, the usual equilibrium version of FDT is obeyed. In 
the low Xo region, the above picture changes, as shown in the bottom panels Fig.^ after 
an early transient the response, A/i, is negative and, thus, X is negative, which would 
asymptotically correspond to a negative "effective temperature" . 

V. CONCLUSIONS 

In conclusion, the present paper has dealt with the understanding of the "jamming" tran- 
sition in granular media via the introduction of models from standard statistical mechanics 
to describe these "non-thermal" systems (the IFLG and Tetris, 

They exhibit a logarithmic compaction when subject to gentle shaking in presence of 
gravity, a compaction extremely close to what is experimentally observed in granular packs 
0]. The presence of such slow dynamics is linked to the existence of a "jamming" transition 
in granular media. In facts, in the high density region, at p^, the models undergo a structural 
arrest where grains self-diffusivity becomes zero 

Self- diffusion suppression at pm signals that, above such a density, it is impossible to 
obtain a macroscopic rearrangement of grains without increasing the system volume, a fact 
interestingly similar to the phenomenon of the Reynolds dilatancy transition P]. 

A very important fact is that results from the present models are in excellent agreement 
with known experimental ones, and their many new outcomes must be still experimentally 
investigated. We have also discussed the intriguing connections of granular media with others 
materials, as glassy systems or spin glasses, where geometrical disorder and frustration play 
a crucial role. 
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FIGURES 






FIG. 1. A schematic picture of the two kind of frustrated lattice gas models described in the 
text. Left: the Tetris model. Right: the Ising Frustrated Lattice Gas, IFLG. Straight and dashed 
lines represent the two kind of interactions eij = ±1. Filled circles are present particles with 
"orientation" Si = ±1 (black/white). 



13 




FIG. 2. The compaction of granular media at low shaking amplitudes. Experimental data 
from Knight et al. (squares) and MC data about IFLG (circles) rescaled according eq. (|^. 
Insert: density p(r, xo;tn) from MC data as a function of tap number t„, for tap amplitudes 
xq = 0.001,0.01,0.05,0.1 (from bottom to top) and duration r = 3.28 • 10^. The superimposed 
curves are logarithmic fits from eq. (^). 



14 




r 

FIG. 3. The static bulk density, p{T), of the IFLG model as a function of the vibration ampli- 
tude, r, in cyclic vibration sequences. The system is shaken with an amplitude F which at first 
is increased (filled circles), then is decreased (empty circles) and, finally, increased again (filled 
squares) with a given "annealing-cooling" velocity 7 = AF/r. Here we fixed 7 = 1.25 10"^. The 
upper part of the cycle is approximately "reversible" (i.e., empty circles and filled squares fall 
roughly on the same curve) . The data compare rather well with the experimental data of Novak 
et al.. F* is approximately the point where the "irreversible" and the "reversible" branches meet. 
Tg signals the location of a "jamming transition". 
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FIG. 4. Main frame: As in Fig.^, the density, p(T), is plotted as a function of the vibration 
amphtude, F for three values of the "cooling" velocity 7. Here only the descending parts of 
the cycles are shown. As in glasses, a too fast cooling drives the system out of equilibrium. 
The position of the shoulder, Tg{'y), schematically individuates a "jamming transition". Inset: 
numerical estimate, in the IFLG model of the dependence of Fg(7) (circles) and F*(7) (squares) on 
the cooling rate 7. Superimposed are logarithmic fits in analogy to those for the glass transition 
temperature, ^^(7), in glasses [4]. When 7 ^ we roughly find Fg(0) = F*(0). 
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FIG. 5. The two time density-density correlation function, {C{t,t') — Coo)/(l — Cqo), as a 
function of the scahng variable a = ln[(t + ts)/r]/ln[(t' + t<j)/r]. Scaled on the same master 
function are data from both models considered in the present paper (Tetris, squares, and IFLG, 
circles) for F = — l/ln(xo) with xq G [10"'^, 10^"*^]. Inset: The correlation C{t,t') for the Tetris at 
F = 0.22 (or xo = 0.01) as a function of t - t' for t' = 10^, 2 10^, 10^, 10"^. 
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FIG. 6. The figures on the left show, as a function of t — t^, the average height difference, 
Ah{t,tw) = hi{t,tw) — ho{t), of a reference system shaken at a given xq and a rephca perturbed 
after by shaking at xq + Axq (Axq = .002). On the right, a check of the generahzed fluctuation 
dissipation relation of eq. (|5|) is given. The integrated response, Ah{t,tw), is plotted as a function 
of the displacement of the reference system, B(t, t^). Systems are shaken at different "amplitudes" 
^0 i^o = 0.8 top, Xq = 0.5 bottom), with replicas perturbed after a = 370. At low xq, negative 
responses appear, but in agreement with eq. (|5|), Ah is asymptotically still approximately linear in 
B. 
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